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Using the 3-telescope IOTA interferometer on Mt. Hopkins, we report results 
from the first near-infrared (A = 1.65/zm) closure-phase survey of Young Stellar 
■ Objects (YSOs). These closure phases allow us to unambiguously detect depar- 

tures from centrosymmetry (i.e., skew) in the emission pattern from YSO disks 
on the scale of ~4 milliarcseconds, expected from generic "flared disk" models. 
Six of fourteen targets showed small, yet statistically-significant, non-zero closure 
phases, with largest values from the young binary system MWC 361-A and the 
(pre-main sequence?) Be star HD 45677. Our observations are quite sensitive 



O 



to the vertical structure of the inner disk and we confront the predictions of the 
"puffed-up inner wall" models of Dullemond, Dominik, and Natta (DDN). Our 
data support disks models with curved inner rims because the expected emission 
appear symmetrically-distributed around the star over a wide range of inclination 
angles. In contrast, our results are incompatible with the models possessing verti- 
ca/inner walls because they predict extreme skewness (i.e., large closure phases) 
from the near-IR disk emission that is not seen in our data. In addition, we also 
present the discovery of mysterious H-band "halos" (~5-10% of light on scales 
0.01-0.50") around a few objects, a preliminary "parametric imaging" study for 
HD 45677, and the first astro metric orbit for the young binary MWC 361-A. 
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Subject headings: accretion disks - instrumentation: interferometers - tech- 
niques: interferometric - stars: binaries: spectroscopic - stars: pre-main se- 
quence - stars: individual (AB Aur, HD 45677, HD 144432, MWC 166, MWC 
275, MWC 297, MWC 342, MWC 361, MWC 480, MWC 614, MWC 863, MWC 
1080, RY Tau, vl295 Aql) 

1. Introduction 

The study of accretion disk evolution in Young Stellar Objects (YSOs) promises to 
directly link the fields of star formation with planet formation. Spatially-resolved, multi- 
wavelength observations can directly probe the dust and gas distribution on solar system 
scales, allowing precise testing of disk models while simultaneously providing the crucial 
physical parameters for planet formation initial conditions. Young Jupiter-mass planets 
embedded in natal disks should even be detectable in the future. 

Already observations from long-baseline interferometers (IOTA, PTI, Keck) have empir- 
ically determined the near-infrared (NIR) emission sizes of YSOs, (e.g., Millan-Gabet et al. 
1999, 2001; Akeson et al. 2000; Eisner et al. 2003), finding the sizes to be strongly linked 
to the luminosity of the central stars (Monnier & Millan-Gabet 2002; Monnier et al. 2005). 
These results are qualitatively explained by a model where the NIR emission arises entirely 
from the hottest dust at the inner edge of an accretion disk, and the measured sizes corre- 
spond to the dust-destruction radius around the young stars (Tuthill et al. 2001; Natta et al. 
2001). Recently, Eisner et al. (2004) reported convincing signs that NIR dust emission from 
many intermediate-mass Herbig Ae/Be stars are elongated, implicating disk-like structure 
as expected from star and planet formation theory. 

The simple detection of "elongation" by an interferometer, however, does not well- 
constrain the geometry of the emission region. For instance, current data can not distinguish 
between a) a centrosymmetric emission pattern as expected from a perfectly flat, inclined 
disk, b) elongated emission arising preferentially from one side of the disk (as expected from 
optical depth effects associated with flaring; Malbet et al. 2001), or c) hybrid disk+halo mod- 
els with multiple emission components (Vinkovic et al. 2006). All current results are based 
on measuring interferometer fringe amplitudes, while interferometric fringe phase informa- 
tion is required to unambiguously detect deviations from simple symmetries. Ultimately, the 
phases and amplitudes are all needed to permit model-independent imaging using aperture 
synthesis techniques (Thompson et al. 2001). 

While atmospheric turbulence corrupts the direct measurement of fringe phase, inter- 
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ferometrists using three or more telescopes can measure the "closure phase" (CP), a phase 
quantity that is unspoiled by telescope-specific phase errors (Jennison 1958; Monnier 2000). 
Recently, the Infrared Optical Telescope Array (IOTA) became the third facility to suc- 
cessfully achieve closure phase operation (Traub et al. 2003; Monnier et al. 2004) and is 
the only such facility that has demonstrated the capability to study Young Stellar Objects 
(Millan-Gabet et al. 2006b). 

Here, we expand the original Herbig Ae/Be survey work of Millan-Gabet et al. (2001) 
using three IOTA telescopes simultaneously, tripling the rate at which visibility measure- 
ments are collected and allowing closure phases to be measured for the first time for a 
sample of YSOs (see Millan-Gabet et al. 2006b, for the first results for the prototype Her- 
big Ae star AB Aur). This work has been made possible by the advanced beam combiner 
IONIC3 (Berger et al. 2003) which exploits integrated optics technology developed for the 
telecommunication industry. 

This article presents a large volume of data and analysis and is organized in the following 
sections. First, we describe the interferometer observations and present the data results in 
a summary form. Second, we analyze the spectral energy distribution and visibility data 
together for each target in order to estimate the NIR emission size; this step is crucial for 
interpreting the closure phase results. Sections 4 & 5 contain the most novel aspects of this 
paper. Here we outline our method for using closure phases to probe the amount of skew in 
the YSO disk emission. Then, we apply our method to quantify the skewed disk emission 
and discuss the implications on current disk models. Lastly, we treat the special case of 
HD 45677 and report the first astrometric orbit for the YSO binary MWC 361-A. 

2. Observations and Data Reduction 

2.1. Description of IOTA Observations 

All the data presented herein were obtained using the IOTA (Infrared-Optical Telescope 
Array) interferometer (Traub et al. 2003), a long baseline interferometer which observes at 
visible and near-infrared wavelengths. Located on Mt. Hopkins (Arizona), IOTA is operated 
by a consortium of institutions, most notably the Smithsonian Astrophysical Observatory 
and the University of Massachusetts at Amherst. The three 0.45-m telescopes are movable 
among 17 stations along 2 orthogonal linear arms (telescopes A & C can move along the 
35-m northeastern arm, while Telescope B moves along the 15-m southeastern arm). By 
observing a target in many different array configurations, IOTA can synthesize an aperture 
35mxl5m (corresponding to an angular resolution of ~5xl2 milliarcseconds at 1.65yum). 
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First results using closure phases have been published by Monnier et al. (2004) and 
Kraus et al. (2005) using the IONIC3 combiner on IOTA. We only just introduce the IONIC3 
combiner here and refer the interested reader to a recent description in Berger et al. (2003); 
an engineering paper with detailed description of the optical component and its performance 
is in preparation (Berger et al. 2006). Light from each telescope is focused into a dedicated 
single-mode fiber and the three fibers are aligned using a silicon v-groove array mated to 
planar waveguides on the integrated optics (IO) device. The optical circuit acts to split 
the light from each fiber /telescope before recombining each telescope pair (AB, BC, AC) at 
three IO couplers. This "pair-wise" combination scheme leads to six interferometer channels 
(two for each baseline) and the fringes are detected using a sensitive HgCdTe (Rockwell 
PICNIC) detector (Pedretti et al. 2004). Varying coupling efficiencies into the fibers (due 
to seeing variations and tip-tilt control errors) cause the system visibilities to vary in time, 
however this effect can be directly calibrated using measurements of the IONIC3 (optical) 
flux transfer matrix which relates the six output channels to the three input channels for 
incoherent light (explained further in Monnier et al. 2004). 

The interference fringes are temporally-modulated on the detector by scanning piezo 
mirrors placed in two of the three beams of the interferometer. A typical single obser- 
vation consists of 200 scans obtained in ~4 min, followed by calibration measurements of 
the background and single-telescope fluxes (important for characterizing the IO beamsplit- 
ters/couplers). Target observations are interleaved with an identical "calibration" sequence 
obtained on an unresolved or partially-resolved star, serving to calibrate the interferometer's 
instrumental response and effects of atmospheric seeing on the visibility amplitudes. The 
target and calibrator sources are typically separated on the sky by 5-10 degrees and are ob- 
served 10-20 minutes apart; these conditions ensure that the calibrator observations provide 
a good estimate of the instrument's transfer function. 

The target sample was selected from the catalog of Herbig Ae/Be members by The et al. 
(1994). Eleven of the 14 targets come from Table 1 of bona fide Herbig Ae/Be members or 
candidate members. HD 45677 is contained in Table 3 "Extreme emission line objects" and 
MWC 342 is found in Table 4a "Other early type emission line stars with IR excess" ; both 
of these objects are discussed individually in §3.2. Lastly, we have included one object not 
listed in the The catalog, RY Tau. RY Tau (F8) is often considered too late-type to be a 
Herbig Ae/Be star and too early to be a T Tauri star. However, the inner disks of T Tauri 
and Herbig Ae stars seem to share common physics (Millan-Gabet et al. 2005; Monnier & 
Millan-Gabet 2002) and we so have included this target in our sample. We note that the 
IOTA Interferometer with the IONIC combiner has a limiting H-band magnitude of ~7, and 
this limitation played a major role in our selections. 
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The science targets, along with their fundamental parameters, are listed in Table 1. All 
the science targets presented here were observed using a standard H band filter (Ao = 1.65/im, 
AA = 0.30/im). Full observing information for our science targets can be found in Tables 2 & 
3, including dates of observation, interferometer configurations, and calibrator details; note 
that the calibrators are all nearly unresolved on our baselines and the size uncertainties do 
not generally significantly contribute to the calibration uncertainties. Figure 1 shows the 
(u,v)-coverage obtained during each epoch for each science target. In the next section, we 
discuss our analysis of these data and the techniques used to process raw data into calibrated 
squared- visibilities (V 2 ) and closure phases (CP). 



2.2. Data Reduction 

Reduction of the IONIC3 visibility data was carried out using custom software similar 
in its main principles to that described by Coude du Foresto et al. (1997, for the FLUOR 
experiment), developed using the Interactive Data Language (IDL). In short, we measure the 
power spectrum of each interferogram (proportional to the target V 2 ), after correcting for 
intensity fluctuations and subtracting out bias terms from read noise, residual intensity fluc- 
tuations, and photon noise (e.g., Perrin 2003). We require that interferograms are detected 
for at least two of the three baselines in order to assure a good closure phase measurement; 
this condition is almost always met thanks to a realtime fringe packet "tracker" (discussed 
further below). Lastly, the data pipeline applies a correction for the variable flux ratios for 
each baseline by using a flux transfer matrix (e.g., Coude du Foresto et al. 1997; Monnier 
2001). We have studied our calibration accuracy by extensive observations of the binary star 
A Vir (Zhao et al. 2006). For bright stars (H mag <5), we have validated 2% calibration er- 
ror for V 2 (corresponding to 1% error in visibility). However, most YSOs are faint requiring 
slow scan rates and exposing our calibration to seeing changes. Under most conditions for 
faint stars, we find a calibration error of ~5% in visibility-squared (e.g., 2.5% in visibility) 
although errors of £10% are occasionally encountered during poor seeing conditions. For 
very resolved objects, the calibration error is dominated by imperfect subtraction of the 
power spectrum background resulting in a AV 2 = 0.02 noise floor. To be conservative, we 
adopted a uniform 5% calibration error (with AV 2 = 0.02 noise floor) for all V 2 in this paper, 
although this overestimates the error for many epochs (Millan-Gabet et al. 2006b, adopted 
similar values for AB Aur). 

In order to measure the closure phase, a fringe-tracking algorithm was applied in real- 
time while recording interferograms (Pedretti et al. 2005), ensuring that interference occurs 
nearly simultaneously for all baselines. We followed the method of Baldwin et al. (1996) for 
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calculating the complex triple amplitude in deriving the closure phase, explicitly guarantee- 
ing that the fringe frequencies for each triple product also "close" [vab + v bc + v ca = 0). 
One calibration step of Baldwin et al. (1996) we do not need to apply is a photon-noise bias 
correction to the triple product (e.g., Wirnitzer 1985), which is only necessary for "all-in-one" 
combiners. 

Pair-wise combiners (such as IONIC3) can have a large instrumental offset for the closure 
phase which is generally calibrated by using a point-source calibrator. The instrumental 
closure phase of IONIC3 drifts less than 1 degree over many hours (a remarkable fact given 
the unstabilized IOTA environment whose temperature can drift by >10 K during a night), 
owing to the miniature dimensions of the 10 component. Also, chromaticity effects limit 
our absolute precision when the calibrator and target are not of the same spectral type due 
to different effective wavelengths when using the broadband H filter. Recent engineering 
tests indicate that the instrumental closure phase (<E>op) of the current IONIC3 combiner 
varies little (^0.5°) between between a hot star (B8) and a cool star (M3) when using the 
broadband H filter. We minimize these errors through repeated calibration of our targets 
with calibrators of similar spectral type. In light of this, we have adopted a systematic error 
of A$c"p = 0.5° for the measurements presented here to represent the residual correction for 
wavelength dependence. This error is not applied in our analysis (it is not random and should 
be fixed for each target/calibrator pair), but rather establishes a minimum closure phase 
magnitude in order to be confident of detecting intrinsic non-zero closure phase (i.e, skew) in 
a target. As a cross-check we note that our targets are generally redder than our calibrators, 
yet we have found no systematic offset in our measured YSO closure phases, consistent with 
expectations from the engineering study. Lastly, we established the sign of our closure phase 
(baseline triangle connecting telescopes A— >B^C^A) using test measurements of Capella 
and Matar in comparison to published orbits (Hummel et al. 1994, 1998). 

In Figure 2, we present all the V 2 data from this survey work, where each epoch of 
observation is denoted by a unique plot symbol. For this and subsequent data analysis, 
the V 2 data from a given configuration has been averaged in the uv-plane (on the scale of 
uv-distance 4-meters) to improve signal-to-noise. We have compared our V 2 results with 
those from previous workers and found good agreement within quoted uncertainties. Note 
that we report here the first ever infrared interferometry results for the targets MWC 342 
and HD 45677. 

We present all the closure phase results in Figures 3 and 4, where the data has been 
plotted against the hour angle of the observation and the array configuration is noted in 
the legend. As was done for the V 2 data, the closure phase data was averaged in the uv- 
(hyper)plane before plotting. 
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All V 2 and closure phase data are available both in averaged and un-averaged form 
from the authors 1 ; all data products are stored in the FITS-based, optical interferometry 
data exchange format (OI-FITS), recently described in Pauls et al. (2005). 

3. Analysis of Visibility Data 

3.1. SED Decomposition and Size Estimation 

In this section, we carry-out the (now) standard procedure of YSO size estimation, first 
outlined by Millan-Gabet et al. (2001). While most of targets have already been observed 
by either IOTA, PTI, or Keck, we have much better (u,v)- coverage than previous studies as 
well as a few new targets (MWC 342, HD 45677), necessitating a fresh visibility and disk 
size analysis. 

Visibilities from YSOs are traditionally interpreted by fitting simple emission models, 
such as Gaussians and rings (e.g., Millan-Gabet et al. 2001; Monnier & Millan-Gabet 2002), 
after correcting for the unresolved stellar emission. For each target, we have collected recent 
visible and near-infrared photometry from the literature to create a spectral energy distribu- 
tion SED (see Table 1 for list of specific photometry references for each object). Using the 
appropriate Kurucz model (Kurucz 1979) for each star (from literature spectral type), we fit 
the SED as a combination of the star plus dust, the latter being simply modelled as a single- 
temperature blackbody; the reddening is another free parameter of this model important to 
match the blue/visible portion of the spectrum. After fitting, we can estimate the fraction 
of flux at 1.65yum arising from the dust component compared to the star. Because YSOs are 
photometrically variable, we included photometry from multiple epochs to illustrate the level 
of variability and this uncertainty dominates the decomposition error budget. Figures 5 to 
18 show the results of the SED decompositions and our estimates of the dust fraction (with 
errors) are tabulated in Table 4. A more detailed description of this algorithm can be found 
in Monnier et al. (2005). 

Next, the V 2 data are fit using a two-component model consisting of the unresolved 
stellar source (flux fraction constrained from SED fitting) and a thin circularly-symmetric 
ring of NIR emission (Gaussian cross-section, 25% thickness) 2 . In some cases, a fully-resolved 



1 Currently, all data can be found at the data archives section of the Optical Long Baseline Interferometry 
News (OLBIN) website, http://www.olbin.jpl.nasa.gov 

2 While our new observations have better (u,v)-coverage than previous work, the relatively short baselines 
of IOTA coupled with the 5% calibration errors on the visibility-squared forbid us from definitively detecting 
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extended component ( "halo" ) was included if evident from the data. These results are found 
in Table 4 and the right panels of Figures 5-18. Note that this "halo" need not be of the 
kind described by Vinkovic et al. (2006), but could arise form scattering off the surface of 
the disk, or from localized emission within even a few AU of the star (Millan-Gabet et al. 
2006b). 

In order to interpret the presence or absence of "halos," one must understand the field- 
of-view of modern infrared interferometers. The most important limitation on field-of-view 
(FOV) is restricted by the field of an individual telescope. For fiber-optic based systems 
(Keck, PTI, IOTA3-IONIC), this FOV corresponds to the diffraction-limited beam which 
is 1.0" for IOTA3-IONIC at H-band, 1.4" for PTI at K-band, and 0.055" for Keck at K- 
band. Note the original IOTA2 experiment (Millan-Gabet et al. 2001) had a "free-space" 
combiner which had a 2-3" FOV. Thus, light beyond this scale does not make it into the 
beam combiner and is completely neglected. The inter ferometric FOV is usually limited 
by bandwidth-smearing. This means that light farther than A ^ B (B is the interferometer 
baseline, A is the central observing wavelength, and AA is the bandwidth) from the central 
source becomes incoherent and does not produce fringes at the phase center (zero optical 
path length difference). In practice, this is about 6x the fringe spacing, corresponding to 
about 60 mas for IOTA (H-band), 25-30 mas for PTI & Keck (K-band). In this latter case, 
light outside the interferometric FOV is detected but is not coherent and acts to lower the 
observed visibility. See Monnier (2003) for further discussion of modern optical and infrared 
interferometry. 

Before discussing individual objects, we note that evidence for large-scale (0.01-0.50") 
"halos" was found for AB Aur (as found by Millan-Gabet et al. 2006b), MWC 275, MWC 297, 
and RY Tau. Also, we surprisingly found that the NIR disk emissions for HD 45677 and 
MWC 614 were completely resolved by the 38 m IOTA baselines; this allowed a unique 
comparison of the predicted stellar flux contribution (from SED fitting) with the observed 
values as well as search for the tell-tale signature of ring-like emission compared to more 
smooth Gaussian-like emission. We discuss unique aspects of some targets below. 

3.2. Individual objects 

Many objects need special attention and specific discussion. In this section we briefly 
describe unusual and interesting characteristics of targets in our sample. 



large-scale elongations (for most sources) such has been reported using the longer-baseline Palomar Testbed 
Interferometer (e.g., Eisner et al. 2004). 
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RY Tau. This star is variously spectral- typed as Kl (Akeson et al. 2005b) to F8 
(Mora et al. 2001), an intermediate object between the Herbig Ae/Be and Classical T Tauri 
stars. Adopting F8 and fitting to the SED, we find a stellar luminosity of ~5 L Q , a few 
times lower than reported by other workers. The quality of the RY Tau visibility data is 
particularly excellent owing to good seeing, and we take seriously the low visibility data 
point at short baselines seen in Figure 5. In this figure, we present fits both assuming a 
10% extended "halo" emission and with no halo emission (both results are also reported in 
Table 4). Clearly, halo emission fits data much better and gives a ring diameter consistent 
with K band results of Akeson et al. (2005b, ring diameter 2.60±0.14 mas), although this 
high level of H-band extended emission seems incompatible with their reported 2% scattered 
light at 0.9/im. 

AB Aur. This object was discussed in detail in Millan-Gabet et al. (2006b). Here, we 
present an independent size analysis, but use the same IOTA dataset. 

HD 45677. Here we report the first long-baseline interferometry data for this target. To 
aid interpretation, we obtained recent visible photometry from the MDM Observatory in 2004 
November and used NIR photometry from 2MASS. Surprisingly, the disk emission appears 
fully resolved at the longest IOTA baselines, a direct measurement of 62±3% dust fraction, 
in reasonable agreement with the value derived from the SED decomposition (65±15%). 

Interestingly, the visibility data does not show the expected shape for ring-like emission 
(specifically, large "ringing" at long baselines due to the ring thinness). Although we were 
able to find an acceptable fit for a ring model, the apparent constant visibilities at long 
baselines motivated additional fitting of a Gaussian emission model and this result, along 
with the best-fit ring model, is included in Table 4 and Figure 8. As is obvious from this 
figure, additional data at slightly longer baselines can rule out or definitively establish ring- 
like emission (as opposed to smooth, Gaussian emission) for this object. The distinction 
may be critical to decide between envelope and inner-rim models of the NIR emission (e.g., 
Vinkovic et al. 2006). 

Here, we first note that the emission appears quite elongated and we additionally find 
large non-zero closure phases. This target is clearly very different from the others in our 
sample, and a more complex model is developed in §6.1 after a more general discussion of 
closure phases. 

MWC 166. Millan-Gabet et al. (2001) found that there is very little, if any, IR excess 
at H band (12 ± 20%) for this high-luminosity target (B0). However, these workers did find 
the long-baseline visibility to be resolved (V~90%), after correcting for 22% emission from 
nearby companion (also B0 spectral type at separation ~0.6", Corporon 1998). We have 



-10- 



obtained new visible photometry at MDM Observatory which, when combined with 2MASS 
data, also suggest a small IR excess (10% ± 10%). 

Our new visibility data also confirm the slightly resolved visibility at long-baselines 
(V~85%), although this is before correction for the known companion. Because the IONIC3 
combiner uses single- mode fibers (unlike the original IOTA two-telescope survey), the com- 
panion flux is more attenuated compared to the primary star, although the exact amount is 
difficult to predict (depending on the seeing and the optical alignment). At short baselines 
(15 m), the observed visibility is ~90% - this limits the companion's contribution to <10% 
of total H band flux. In order to account for this, we performed fits including both 0% and 
10% contribution from companion, while also constraining the dust fraction to be between 
1-20%. The allowed ring diameters span the range 4-15 mas, with a best-fit for 80% primary 
star, 4% secondary star (incoherent flux), and 16% disk emission using a ring diameter of 
7.5 mas. Some representative fits are shown in Figure 9. 

MWC 863. A binary companion was detected by Corporon (1998), a T Tauri star (K4) 
at 1.1" and AH ~2.24. We see no evidence for this companion in our visibility data, since 
the field-of-view of an IONIC3 single-mode fiber is only ±0.5". Although the companion flux 
was taken into account in our estimate of the dust fraction in SED fitting, the SED fit is 
still poor in the near-IR region. While our measured V 2 is similar to that found by Millan- 
Gabet et al. (2001), our ring diameter estimate is much smaller because we did not apply 
the same visibility "correction" to compensate for the incoherent flux of the companion (our 
result here is consistent with the recent result using the longer baseline Keck Interferometer 
at K band; Monnier et al. 2005). Also, our luminosity determination is significantly lower 
here than recently reported in Monnier et al. (2005), due to selection of more recent optical 
photometry. Better coeval photometry is needed for this target in order to estimate stellar 
and disk properties precisely. 

MWC 297. V 2 data are better fit using a 5% extended "halo" emission, although 
the fit is reasonable assuming no halo (see Figure 13). The best-fit values (and uncertainty 
estimates) in Table 4 include a 2.5±2.5 % estimate for halo emission. This is one of the most 
resolved disks in our survey and, while our uv-coverage is limited, there is slight indication 
of elongation along PA~100°; more data are needed to confirm this. 

MWC 614. The nearly flat visibility curve indicates the dust emission is entirely 
resolved even at short baselines, and that the star makes up 61% of the flux, consistent with 
SED decomposition estimate of 63±10 %. We can convert this to a lower limit on the dust 
envelope size, ring diameter >13 mas (best-fit Gaussian FWHM 15.9 mas; see Figure 14). 
However, as for HD 45677, we see no telltale sign in the V 2 data for a ring-like structure 
which has the distinctive Bessel function shape in Fourier-space (e.g, Gaussian profile is a 
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better fit than ring profile). This target would be ideal source for adaptive optics on an 
8-m class telescope. The large ring size derived here is anomalous for the Herbig Ae size- 
luminosity relation (Monnier et al. 2005). This could be partially-explained if there is a stellar 
companion within about 1" of the primary star, since this would lead to an overestimate of 
the dust ring size in our study. 

MWC 342. We present first long-baseline interferometry data on this target. While our 
data indicates elongated emission, our very limited data set does not allow us to meaningfully 
quantify this. We note that the evolutionary state of MWC 349 is unclear; Miroshnichenko 
et al. (1999a) suggests this star, as well as HD 45677, are possibly X-ray binaries like CI Cam. 

MWC 361-A. A 6" companion to MWC361-A was first discovered by Li et al. (1994). 
While this object does not affect our measurements, the close companion (separation ~18 mas) 
discovered by Millan-Gabet et al. (2001) obviously is easily detected and appears to be re- 
lated to the ~3.7 yr periodic variations in the Ha emission seen by Miroshnichenko et al. 
(1998). Indeed, Pogodin et al. (2004) recently reported a detailed spectroscopic orbit (more 
discussion of this in §6.2). Here in Figure 17, we have plotted the V 2 of primary star after 
correcting for the fainter companion contribution (made possible only after the full orbit 
analysis of §6.2). Note that this ring diameter estimate assumes the companion star is it- 
self unresolved, a reasonable assumption based on the observed size-luminosity relations. 
Interestingly, the NIR ring diameter found here is undersized compared to the normal size- 
luminosity relation for Herbig Ae/late Be stars, as also found for some other early Be Herbigs 
(Monnier & Millan-Gabet 2002; Eisner et al. 2004). 

MWC 1080. There seems to be a fair amount of uncertainty on the spectral type 
of this star, ranging from BO-2 (Cohen & Kuhi 1979) to AO-3 (Yoshida et al. 1992), and 
distances of 1 to 2.5 kpc. Here, we adopt the earlier spectral type and a distance of 2.2 kpc 
(Grankin et al. 1992). In principle, our visibility data should be sensitive to the "wide" 
binary (Corporon 1998) seen with a separation of ~0.8" (AH ~ 2.8), however there is no 
sign of this in our data, presumably due to high brightness ratio and relatively large off-axis 
distance. 



4. Closure Phase as Tool for Studying YSOs 
4.1. Background 

Basic interferometry theory informs us that an interference fringe measured on an in- 
dividual baseline corresponds to a single Fourier component of the brightness distribution 
under scrutiny. Basic Fourier theory further informs us that any brightness distribution can 
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be reconstructed by collection of sufficient Fourier components (i.e., "visibility" data). Thus, 
an "imaging" interferometer, such as the radio Very Large Array (VLA), is characterized 
by a large number of telescopes (corresponding to many simultaneous Fourier components) 
and a range of baselines (for imaging a range of spatial scales). While the VLA accom- 
plishes this using 27 telescopes with a maximum spacing/minimum spacing of £100, current 
infrared interferometers have 3-6 telescopes with maximum/minimum spacings of ^10. Not 
surprisingly, current interferometers have difficulty performing model-independent image re- 
constructions (except for simple scenes such as binary stars). 

The problem is made more vexing for optical interferometry since atmospheric turbu- 
lence scrambles the Fourier phases of the observed interference fringes. In some cases, these 
Fourier (fringe) phases are not important - it can be easily proven that centro symmetric im- 
ages can be completely characterized by the fringe amplitudes alone 3 . From this perspective, 
fringe phases tell us about asymmetric information in an image and, as such, are unneces- 
sary for characterizing simple objects such as (idealized) stellar photospheres or ellipsoidal 
distributions (such as flat inclined disks). 



4.2. What Is the Closure Phase? 

The closure phase method was developed by early radio interferometrists (Jennison 
1958) to calibrate phase drifts in phase-coherent (heterodyne) detection systems, and it was 
first applied to optical data in the context of aperture masking (Haniff et al. 1987). The 
basic idea can be explained by considering a three-telescope array as pictured in Figure 19. 
Consider that turbulent air pockets above aperture 2 will cause equal but opposite effects 
on the two baselines connecting apertures 1-2 and apertures 2-3. Thus, summing up phases 
around a triangle of 3 telescopes cancels out all atmospheric disturbances within the "closed" 
triangle, resulting in a good observable, the closure phase, that is independent of the phase 
fluctuations above the array of telescopes. 

The following derivation mathematically expresses the qualitative ideas of the last para- 
graph. Consider the the closure phase ($^) for a triangle connecting telescopes i, j, and 
k. 

*g* = + $ ? + <^ s (i) 

where, $° bs here is the observed fringe phase for baseline between telescopes i and j. The 



3 By defining the coordinate origin at the center of symmetry, all the Fourier phases become or n 
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observed fringe phase can be seen to be equal to the intrinsic phase ^ (this is what we are 
most interested in) and the difference between the time-varying atmospheric phase delays fa 
associated with individual telescopes. Thus, 

$g> s = (2) 

intrinsic atmospheric 

Putting these last two equations together we can prove the closure phase to be inde- 
pendent of atmospheric phase delays: 



*S = + + (3) 

= + (fa - fa) + + (fa - fa) + + (0; - fc ) (4) 

= + + (5) 

In principle, images can be reconstructed using only visibilities and closure phases (re- 
lated to the technique of "self-calibration" in radio interferometry; Cornwell & Wilkinson 
1981) when additional constraints are imposed (such as limited field-of-view and positivity). 
A more pedagogical description and derivation of the closure phase can be found in Monnier 
(2000) and Monnier (2003). 



4.3. Important Properties of Closure Phases 

The most important property of closure phases follows directly from the properties of the 
Fourier phases already discussed. That is, a centro symmetric (point-symmetric) image will 
always have a closure phase of zero degrees 4 . Thus, detection of non-zero closure phase is a 
robust sign of skewed emission, a definitive deviation from point-symmetry. This remarkable 
and robust property of closure phases is key for studies of young stellar objects, since model- 
independent image reconstructions are today beyond our capabilities (due to limited uv- 
coverage) . 

Another key point is that the interferometer must have sufficient angular resolution to 
resolve the scale over which the emission is skewed: an unresolved source is always point- 



4 Actually, 180 degrees is also possible but not relevant for our case of resolved structure around a dominant 
point-like star. 
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symmetric. Thus, one must match interferometer resolution to the appropriate angular 
scale. 

While the detection of skewed emission is robust, quantifying the degree of skew is very 
model-dependent. When limited uv-coverage forbids direct image reconstruction, one must 
carry-out model fitting. Fortunately, disk theorists routinely produce models for fitting 
to spectral energy distributions and the results of radiative transfer models are suitable 
for comparing to interferometry data (e.g., Malbet & Bertout 1995; D'Alessio et al. 1998; 
Whitney et al. 2003; Lachaume et al. 2003; Akeson et al. 2005a). 

4.4. Skewed NIR disk emission from YSOs 

In this paper, we compare measured closure phases with physically-motivated models 
for the disk emission in order to quantify the level of skewed emission present in our sample. 
While future closure phase surveys with larger telescopes should provide a rich and detailed 
dataset for modeling, we show below that even the limited, first-generation closure phases 
presented here can distinguish between current competing disk models. Here, we contemplate 
two classes of disk models: the DDN disk model and the ad hoc "skewed ring" model. 

44.1. DDN Models 

Dullemond, Dominik, & Natta (2001, DDN) established the current paradigm for Herbig 
disks (and by extension, even new T Tauri disk models; Muzerolle et al. 2003). As initially 
outlined by Natta et al. (2001), the star is surrounded by an optically-thin cavity with a 
gas-only disk. The dusty disk begins when the temperature drops below the dust destruction 
temperature (T ~ 1500K), at which point the midplane becomes optically thick. At the dust 
destruction radius the disk is clearly quite hot and this "puffs-up" the inner rim. The main 
achievement of Dullemond et al. (2001) was to couple the radiative transfer to a (vertical) 
hydrostatic disk structure, finding a physically self-consistent height for the inner rim and the 
outer disk flaring (which match on to traditional flared disk models, e.g. Chiang & Goldreich 
1997). 

In the DDN paradigm, the frontally-illuminated inner wall is perfectly vertical. These 
authors comment that such a situation may not be physically reasonable, and indeed a variety 
of effects may "curve" the inner rim (some of these are discussed in §5.3). In the case of a 
vertical wall, an observer will generally only see the far side of the inner rim wall, since the 
front side will be blocked by outer disk material. This scenario produces maximally-skewed 
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disk emission, boding well for closure phase studies of YSOs. 

Figure 20a shows a synthetic image for the near-IR emission from a DDN disk with 
6 mas inner rim diameter, inclined to our line-of-sight by 30°. The hot inner wall is assumed 
to have uniform surface brightness and we neglect any near-IR emission or scattering from 
the rest of the (cooler) disk. Here we take the half-height of the wall to be ~0.3x inner 
radius (following Dullemond et al. 2001). As expected, the disk emission is highly off-center 
- skewed - with respect to the star. 

We can calculate the expected closure phase signal for this disk model using the IOTA 
interferometer geometry shown in Figure 20b, a typical layout used for the observations in 
this paper. Figure 20c show the predicted closure phase signal for the DDN model (and also 
a rotated version) for A = 1.65/im. Here we also assumed that 65% of the emission came from 
the dust while 35% arises from the star (typical decomposition for our targets). Figure 20d 
shows how the signal is much stronger if the interferometer has 2 x longer baselines - clearly, 
the disk appears more "skewed" with increasing angular resolution. 

44.2. Generic "Skewed Ring" Model 

Later in this paper, we will show that the DDN model predicts large closure phases that 
are not observed. In order to move beyond the vertical wall assumption of the DDN model, 
we will explore the class of "skewed ring" models. Monnier & Millan-Gabet (2002) discussed 
in detail the merits of a ring model for fitting to the near-infrared visibilities of YSOs, and 
this basic model has been used by most current workers in this field (e.g., Millan-Gabet 
et al. 2001; Eisner et al. 2004; Akeson et al. 2005a). The narrow range of temperatures over 
which dust emits at 1.65yum and the empirical results of Tuthill et al. (2001, imaging of 
the LkHalOl disk) both suggest the dust at the inner edge of the YSO accretion disk will 
appear as a thin ring in the near-infrared. Here we adopt a Gaussian cross-section for the 
ring emission, in order to approximate smoother (i.e. curved) inner rim emission. 

Here, we accept this ring model with an additional complication - we modulate the 
ring brightness as a function of azimuth by a sinusoid (essentially, an m = 1 mode). This 
model was recently used for interpretation of the AB Aur disk by Millan-Gabet et al. (2006b). 
Figures 21a & b show two different ring models with a "skew" of 1 and 0.5 (and ring thickness 
of 25% using Gaussian profile). The "skew" refers to the amplitude of the modulation in the 
ring brightness, thus a skew of is centrosymmetric, a skew of 0.5 has a 3-to-l brightness 
contrast between the brightest and faintest (diametrically-opposed) portions of the ring, and 
a skew of 1 has a contrast formally of 00. You will note that the skew 0.5 image is remarkably 
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similar to the image of LkHalOl reported by Tuthill et al. (2001). 

Figures 2 led show the observable closure phases for these models using the same ob- 
serving geometry introduced for Figure 20. The closure phase signal for the skew= 1 model 
is smaller than the DDN model, as expected since the emission is far less skewed. Also, we 
see that the skew= 0.5 model presents even smaller closure phases for the same observing 
conditions. These calculations illustrate clearly how precision closure phase measurements 
of YSOs by today's interferometers can determine the "skewness" of the near-IR emission, 
discriminating between competing disk models (see §5). 

If the skewness in disk emission is a result of a disk viewed at intermediate inclination, 
we might plausibly expect the skewed ring to be elliptical and not circular. This level of 
detail is beyond the current constraints from the interferometer data and would introduce 
another free parameter. Thus, we will not consider this reasonable extension to the "skewed 
ring" model for the majority of our targets (we do introduce such a sophisticated model for 
the special case of HD 45677 in §6.1). 

5. Analysis of Closure Phases 

Except for AB Aur (Millan-Gabet et al. 2006b), none of our targets have previously 
been observed using closure phases and we will concentrate our analysis of this section on 
this observable. We find 8 of the 14 sample targets have closure phases consistent with zero 
within measurement errors. Only AB Aur, HD 45677, MWC 297, MWC 361-A, MWC 614, 
MWC 1080 show statistically-significant non-zero closure phases (at 2-a level). 

In order to interpret this result, we must consider the angular resolution of the IOTA 
interferometer and disk sizes of our targets (measured in §3), since many of our disks are only 
marginally-resolved. Recall from §4 that closure phases encode deviations from centrosym- 
metry in an extended object, closely related to the "skewness" of the emission. Perhaps best 
explained in Lachaume (2003), marginally resolved objects naturally possess small closure 
phases even if the inherent emission is extremely asymmetric. An intuitive explanation of 
this is to consider that any "unresolved" source must be considered to be centrosymmetric. 

Thus, our first priority is to determine whether the finding of small closure phases reflects 
a true lack of skewness in the YSO disk emission, or rather results from from inadequate 
angular resolution. To further this investigation, we have collected the maximum closure 
phase observed for each target in Table 5, along with the disk size measured in §3 in units 
of the fringe spacing (X/B) of the longest IOTA baseline in the appropriate triangle. 



-17- 



We confront the specific predictions of the DDN model in the next section. 

5.1. Confronting DDN Models 

In order to extract quantitative information from the near-zero closure phases of most 
stars in our sample, we will first explore the closure phase predictions of DDN (Dullemond 
et al. 2001) models, described in §4.4.1. 

In order to explore parameter space, we generated closure phase predictions for DDN 
models with a range of inner rim diameters and viewing angles. In each case, we assumed 
the disk made up 65% of the emission at H-band with the remaining 35% arising from the 
(unresolved) star, a typical case for our target sample. We then generated synthetic closure 
phases for the actual observing triangles used by IOTA. The main results of this parameter 
study are shown in Figure 22. 

In Figure 22, we plot the maximum observable IOTA closure phase for DDN disks for 
three different inclinations as a function of the ring diameter. The ring diameter is expressed 
here in units of the fringe spacing X/B of the longest baseline in a given IOTA closure 
triangle. For each skewed ring we sampled all possible position angles on the sky and here 
plot the maximum value of the closure phase magnitude. While we typically do not know 
the sky orientation of a given object, we expect our sample of 14 targets to sample random 
position angles on the sky. 

This procedure quantifies the qualitative results from the DDN examples in Figure 20, 
that well-resolved DDN disks have large closure phase signals. In Figure 22, we can then 
compare the synthetic curves to the maximum closure phases actually observed in our survey. 
To do so only requires us to account for the differing sizes of each target from §3. We exclude 
MWC 361-A from this analysis since the large closure phases for this target likely come from 
the binary companion, not skewness in the disk around the primary star. 

Based on this comparison of DDN model predictions with our closure phase dataset, we 
can draw 2 primary conclusions. 

• The most well-resolved disk sources have smaller closure phases than the DDN model 
calculations, though usually non-zero. This indicates some modest skewness in the disk 
emission but is incompatible with the large skew signal predicted by DDN models. 

• When the ring diameter is less than the ~ | fringe spacing of the longest baseline in a 
triangle, the closure phase signal is strongly suppressed. This applies to approximately 



-18- 



half of our target sample. We will require approximately 2-3 x better angular resolution 
to fully characterize disk skewness in our sample (corresponding to 70- 100m baselines 
at H band or 100- 150m baselines at K band). 

For any given target, the putative disk skewness might be orientated along a position 
angle not sampled by the IOTA interferometer, which is most sensitive to skew in the N-S 
direction (due to non-uniform uv-coverage). Thus, even highly skewed disk emission can 
be "hidden" from the interferometer due to unfavorable observing geometry. However, the 
small observable closure phases for the most resolved disks (MWC 297, MWC 166, MWC 614, 
MWC 480, and RY Tau), when taken as a whole, argue strongly against the vertical wall 
of the inner rim put forward in the DDN model, since this geometry always predicts strong 
skew in disks viewed from all inclination angles. Note that the DDN model (with vertical 
inner wall) also predicts highly elongated disk emission (see Figure 20) even for nearly face-on 
disks, inconsistent with the V 2 results of Eisner et al. (2004). Note that AB Aur uniquely 
appears to show definite skewed emission in its halo (beyond the inner rim), further discussed 
in Millan-Gabet et al. (2006b). 



5.2. Confronting "Skewed Ring" Models 

The analysis of the last section was also carried out using the "skewed ring" model 
described in §4.4.2. The results of this study, along with the observed closure phases, are 
plotted in Figure 23. While skew= 1 disks do generally overpredict the closure phase magni- 
tude, our data can not clearly distinguish between the models with varying "skew" factors. 
Although ad hoc, the skewed ring models are physically motivated and can currently explain 
the observed visibility and closure phase data in our sample, including the the imaging data 
for LkHa 101 and the closure phase data of HD 45677 (taking ellipticity into account; see 
§6.1). While it is critical to develop physical models for the disk emission, the introduction 
of the skewed ring model offers a new and useful tool for describing the emission independent 
of the underlying physical models. 

For some disks, it should be possible to do more detailed model fitting by including 
inclination and position angle information known from CO data or scattered light images 
(e.g., AB Aur, MWC 275). In these cases, the IOTA closure phases can play a more definitive 
role in testing specific models. However, in general, longer baseline closure phases and more 
uniform Fourier coverage will be needed to refine future models of the inner rim geometry. 
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5.3. Discussion 

5.3.1. Implications on inner rim geometry 

The most significant result from this paper is the report of few skewed disks around 
Herbig Ae/Be stars. This is robustly shown by the uniformly small closure phases found 
in our sample of YSOs. When the sample is viewed as a whole, the small closure phases 
have been shown to be incompatible with generic predictions of the original DDN disk mod- 
els; these models possess vertical walls at the inner edge of the accretion disk introducing 
strongly-skewed emission at all viewing angles. Note that this result does not invalidate the 
overall success and applicability of the DDN model which mostly focuses on other aspects of 
the disk structure. While we have shown that "skewed ring" models are better suited for ex- 
plaining the small observed closure phases, we also demonstrate that detailed determinations 
of "skew" required longer baselines than possible with IOTA. 

In the context of today's disk models, our data lend support to inner disk models like 
that of Isella & Natta (2005) which includes a curved inner rim. In this case, the inner 
rim curves away from the midplane due to the expected pressure-dependence of the dust 
evaporation temperature. We plan to explore quantitatively the expected closure phase 
signals of this model with future radiative transfer work, but this is beyond the scope of this 
paper. 

Note that there are other mechanisms that could cause the inner rim to "curve." One 
new idea we suggest is a consequence of dust settling and growth which leads to larger dust 
grains in the midplane than in the upper layers of the disk. Since large grains better radiate 
in the infrared, they are able to exist closer to the star. Smaller grains in the upper layers 
have difficulty cooling and thus evaporate more easily. Thus, the inner rim can be highly 
curved if the average dust size varies significantly as a function of disk scale height in the 
inner rim. This process can cause a much greater curvature in the inner rim compared to 
the effects of gas pressure on the evaporation temperature. Indeed, other recent work have 
found evidence for dust settling and growth in YSO disks (e.g., Rettig et al. 2005; Duchene 
et al. 2004), and the implications on the inner rim structure have only just begun to be 
explored (e.g., Tannirkulam et al. 2005). 

Understanding the geometry of the inner disk is not just important for understanding 
the near-IR emission, van Boekel et al. (2005) showed that 20-30% of the mid-IR emission 
can come directly from the inner rim. Given the intriguing changes in the silicate feature 
observed using the MIDI instrument on VLTI (van Boekel et al. 2004), it is critical to 
understand how the dust properties of the rim might be different than in the surface layers 
of the disk. In addition, high spectral resolution CO observations of near-IR fundamental 
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and first-overtone lines (Najita et al. 2003; Brittain et al. 2003; Blake & Boogert 2004) require 
temperature and density modeling of the inner AU of Herbig disks. Currently, these models 
have not incorporated the puffed inner wall into their calculations, an improvement urgently 
needed in order to correctly interpret the potentially powerful kinematic and temperature 
information available. 



5.3.2. The "Halo" Phenomenon 

By using three telescopes, our survey has measured visibilities on a larger range of 
scales than previous work using only two telescope interferometers. We find evidence for 
intermediate- to large-scale "halos" around AB Aur, MWC 275, MWC 297, and RY Tau. 
For our purposes here, we define halo emission to be any extended emission on scales between 
10-1000 milli-arcseconds making up 5-20% of the near-IR flux. Similar extended halos have 
been reported previously (albeit on 0.5-2.0 arcsecond scales) around some T Tauris and 
Herbigs using infrared speckle interferometry (Leinert et al. 1993, 2001) and even on smaller 
scales probed by lunar occultation (e.g., case of DG Tau; Leinert et al. 1991; Chen et al. 
1992). Similarly-sized extended emission has recently been reported around a sample of 
FU Ori objects using the Keck Interferometer (Millan-Gabet et al. 2006a). 

Existence of these halos is not well-understood theoretically, but could arise from diffuse 
scattered light at large (5-100 AU) radii (as recently discussed by Akeson et al. 2005b). from 
the central star or even from localized thermal emission within a few AU. The origin of the 
halo material is unclear, but plausibly could be from an infalling remnant envelope, dust 
entrained in the stellar wind/outflow, or (in the case of FU Ori) ejected gas & dust from an 
earlier outburst. The connection to the hybrid halo+disk model of Vinkovic et al. (2006) is 
ambiguous and our findings should motivate further work on such multi-component models. 
In the special case of AB Aur (Millan-Gabet et al. 2006b), we detected non-zero closure 
phases from the "halo" component, leading us to conclude the source of the emission was 
localized and within 1-4 AU of the star. For the other "halos" reported in this paper, we 
can only roughly constrain the location of the emission (10-500 milliarcseconds) ; further 
observations with short baselines will be required to elucidate the nature of "halos" around 
Herbig Ae/Be stars. 

Recently, Baines et al. (2006) report interesting "spectro-astrometry" observations of 
many Herbig stars finding evidence for possible binaries around many targets (including 
AB Aur, HD 45677, vl295 Aql, MWC 361-A of this sample). While we caution that a 
few of our "halo" detections may turn out to be due to very close binary companions (10- 
500 mas), most stars in our sample have been surveyed for companions using adaptive optics 
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down to about 100 mas. (e.g., AB Aur, MWC 297, vl295 Aql; Eisner et al. 2004). Similarly, 
the spectro-astrometry data can misinterpret halo emission as evidence for binarity if the 
extended emission itself is skewed (as has been shown for AB Aur; Millan-Gabet et al. 
2006b). Note that scattered disk/halo light can only affect spectro-astrometry results if 
it has a substantially different Ha spectrum than what we see from the star. This could 
happen perhaps in two ways. The Ha emission from the star is probably not uniformly 
distributed, but more concentrated near polar regions and will also have a different dynamical 
signature/line profile - thus, the observer will see different Ha than the disk midplane which 
then scatters the light into our line-of-sight. Second, stellar light must traverse long paths 
just above the disk midplane before scattering into our line of sight, thus there could be 
significant Ha absorption (depending on ionization state). These effects deserve more study 
before settling on an interpretation of the spectro-astrometry data. 



6. Special Treatment for HD 45677 and MWC 361 

6.1. Parametric Imaging of the HD 45677 Disk 

One exception to the general pattern of small closure phases is HD 45677, which is both 
heavily resolved, elongated, and shows strong closure phases. We already fit a symmetric 
ring and Gaussian to the model in §3 and Figure 8, and noted the large fit residuals due to 
elongated and skewed structures evident in the raw V 2 and closure phase data. 

First we fit the V 2 data with an elongated ring model, and the results of this fit are 
shown in the top panels of Figure 24. The direction of elongation (PA 70°) is similar to the 
observed (visible) polarization angle of 60-80° (Coyne & Vrba 1976). However, as expected, 
this model can not explain the relatively large observed closure phases that must arise from 
skewed emission. The results of this fit are tabulated in Table 6. 

In order to explore this, we allowed the elongated ring to be skewed, following the 
same method as for the skewed ring model in §4.4.2. That is, the elongated ring emission 
is modulated as a function of azimuth. The best fit parameters are found in Table 6 and 
the synthetic image is shown in Figure 24 (bottom panels). Interestingly, the skew angle 
(position angle from central star with greatest skewed brightness) is ~90 ° different than 
the major axis of the elongated ring, a geometry that might be expected for a flared disk 
(e.g., Malbet et al. 2001). Given our parametric image and the direction of the skew, we 
would predict that scattering by dust would produce a net linear polarization of PA ~70°, in 
agreement with observations. The x 2 At to the V 2 is poor at the shortest baselines suggesting 
additional complexity missing from our model. 
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The skewed ring model presented here can be viewed as the first attempt to do "paramet- 
ric" imaging of YSO disks using long baseline interferometry. While the number of degrees 
of freedom in our model is still relatively small, we are reconstructing fairly complicated 
image characteristics. Similar quality "images" should be possible for the rest of our target 
sample in the coming years. 

Lastly, we want to comment that this ring image is qualitatively similar to the expected 
results outlined in previous papers (Coyne & Vrba 1976; Oudmaijer & Drew 1999). Based 
on the relatively cool infrared excess and the linear polarization angle, a ring-like, elongated 
morphology was predicted and the position angles found here are entirely consistent with 
the polarization signatures. Combining infrared polarimetry with interferometry could reveal 
unique aspects of the dust in this system (e.g., Ireland et al. 2005). 

We believe that our data support the idea that HD 45677 is not a bona fide Young Stellar 
Object, but rather a more evolved system akin to CI Cam (e.g., Miroshnichenko et al. 2002). 
The disk size is too large to be a typical YSO system (unless we are observing a unique 
transition object). Perhaps a binary interaction in the recent past created the disk-like dust 
distribution around this source (e.g., Miroshnichenko et al. 2005, and references therein). 
Studies of this dust shell may shed light on the mysterious process of dust production around 
some main-sequence Be stars. 

6.2. Astrometric Orbit for the Young Binary MWC 361-A 

The first definitive detection of the MWC 361-A close binary was reported by Millan- 
Gabet et al. (2001) using interferometry. As previously already mentioned, spectroscopic 
evidence has also been accumulating and Pogodin et al. (2004) recently reported a prelimi- 
nary single-lined radial velocity (RV) orbit of the primary with a period of 1341 days. 

There are few simultaneous radial velocity and astrometric orbits for young star binaries 
(e.g., HD 98800-B; Boden et al. 2005) 5 These data are particular valuable because they 
directly constrain the stellar masses and sometimes the system distance, invaluable "ground 
truth" tests for theories of stellar interiors and evolution. 

We report five new epochs of V 2 and closure phase data for MWC 361-A in this paper. 
All these epochs have been fitted for a binary system, with the additional complication that 
the primary is clearly partially resolved. For the treatment in this paper, we have assumed the 
fainter secondary star is completely unresolved as expected for the size-luminosity relations 



5 Eclipsing binary systems are also powerful tools for determining stellar masses (e.g., Stassun et al. 2004). 
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of Herbig stars. We also have independently re-fit the two epochs of data from Millan-Gabet 
et al. (2001) in order to extend our time coverage over ~5 years, more than 1.5 orbital 
periods. 

Table 7 shows our best-fit positions for the secondary component relative to the primary 
for each of the seven epochs. Furthermore, our fitting found a brightness ratio of 6.5±0.5 
(primary/secondary) at H band and a primary uniform disk diameter of 3.6±0.5 mas (we 
analyzed the primary disk size in §3 and Figure 17). There was some evidence for variations 
in the flux ratio with time, as expected for young stars and previously observed for the 
Z CMa system (e.g., Millan-Gabet & Monnier 2002). For about half the epochs (1998 
June, 2003 November, 2004 June, 2005 June) we were able to uniquely constrain all the free 
parameters. For the remaining, we had to fix the brightness ratio and size of the components 
to yield a good positions (due to more limited datasets). In a few cases, we had to use a 
priori information concerning the likely binary position to break position degeneracies. For 
instance, the 1998 September data comes only from one night and multiple solutions exist 
- however, only one solution is sufficiently close to the unambiguous result of 1998 June to 
be considered plausible and consistent with the other epochs. 

Figure 25 shows error ellipses for the position of the secondary relative to the primary for 
all epochs. The position predictions for our best-fit orbit are also shown and the agreement is 
excellent. The orbital parameters (and conservative error estimates) are contained in Table 8 
and compared to the orbital elements of Pogodin et al. (2004). All orbital elements agree 
within observational uncertainties. A future paper will fit the V 2 and closure phase data 
directly (along with radial velocity data) using the orbital elements (as recently done using 
IOTA data for A Vir by Zhao et al. 2006). This analysis will yield a joint orbital solution 
with significantly- smaller parameter errors, however this analysis is beyond the scope of this 
paper. 

With the combined RV and astrometric orbital solution, we can determine physical 
characteristics of the MWC 361- A system. We had the Hipparcos intermediate data re- 
analyzed using the new orbital parameters in an attempt to improve the parallax and 
constrain the photocenter orbit (in the method of Pourbaix & Jorissen 2000). The new 
result n =2.76±0.68 mas (d= 360lyQ° pc) varies little from the original catalog result 
(tt =2.33±0.62 mas;d= 430+go° pc). Using the orbital period (3.74±0.06 yrs; weighted- 
average between the independent RV and astrometric determinations), and (angular) semi- 
major axis (15.14±0.70 mas), we can calculate the total system mass M 1 + M 2 = 10.4; 20 9 5 
using the original Hipparcos measurement (these errors represent 1-a confidence intervals 
using a Monte Carlo method for error propagation). Unfortunately, while these preliminary 
masses are consistent with our expectations for a binary system of early B stars, the large 
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distance errors propagate into large mass uncertainties, rendering these estimates useless for 
constraining any models. It will be crucial to extract a double-lined radial velocity orbit in 
order to determine an orbital parallax, resulting in accurate and precise mass estimates for 
both components. 

We note that Pogodin et al. (2004) measured the mass function f(M 2 ) for the secondary 
from the single RV orbit. This result (f(M 2 ) = 0.175 ± 0.035) was interpreted using an 
assumed inclination i = 70°, consistent with our measurement of i = 65 ± 8°. Thus, we 
can concur with their analysis that MWC 361-A appears to be a close binary consisting of 
an early Be star, likely hosting the resolved disk of §3, and a lower- mass, late Be star as 
secondary. 

The periastron distance of the MWC 361-A components is only ~4 AU (due to somewhat 
eccentric e ~ 0.3 orbit), thus the accretion disk around the MWC 361-A primary should be 
severely truncated (perhaps this relates to recent difficulty in fitting the far-infrared SED, 
Elia et al. 2004). Certainly, MWC 361-A provides a fascinating laboratory for studying 
binary star formation and is suitable for future studies by longer baseline interferometers 
that can both measure orbital parameters at high accuracy but also image the accretion 
disks around both stellar components. It is extremely rare to find young stars of such early 
spectral types in a binary system, and MWC 361-A promises to be a new anchor point for 
models of stellar interiors and evolution, as well as for disk evolution in binary systems. 



7. Conclusions 

Here we have presented the results of the first closure phase survey of Young Stellar Ob- 
jects using an infrared interferometer. We find most YSO targets show no signs of skewness 
in their emission on the scale of ~4 milliarcseconds. Our data is incompatible with vertical 
wall models of the puffed-up inner rim (e.g., Dullemond et al. 2001). A "skewed ring" model 
was developed in order to quantify our results and the observed small closure phases favor 
inner dust rims which are smoothly curved (e.g., Isella & Natta 2005). At least a crude 
understanding of the inner AU of YSO disks is crucial for interpreting a host of other YSO 
observations, including spectrally-resolved CO data, spectro-astrometry, mid-IR sizes and 
spatially-resolved silicate feature minerology. Baselines significantly longer than the maxi- 
mum IOTA length of 38 m (i.e., from CHARA and VLTI) will be required to investigate this 
further for specific sources. 

An independent and complete reanalysis of new H band V 2 data found fresh evidence 
for extended "halos" of emission with about 5-10% of the emission for a few targets (see 
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also Leinert et al. 2001). Halos appear to be quite common in many young disk systems 
(T Tauri, Herbig, FU Ori), although the physical origin of this large-scale emission appears, 
as of yet, not understood. These new data may re-open the debate concerning the relative 
importance of disks, halos, envelopes, localized thermal emission, and the applicability of 
hybrid models. 

We also took the first steps towards imaging young stellar objects by studying the 
unusual system HD 45677, finding this Be star to be surrounded by an elongated and highly- 
skewed dust ring. We agree with earlier suggestions that HD 45677 may not be a young 
object, but rather a member of the enigmatic class of main sequence hot stars with large IR 
excess. 

Lastly, we calculated the first astrometric orbit for the short-period binary MWC 361- 
A, finding good agreement with recent spectroscopic work. We found mass estimates in 
line with expectations and discuss the exciting science potential from future studies of this 
unique system. 
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Table 1. Basic Properties of Targets 



Target 


RA (J2UUUJ 


Dec (J2UUUJ 


V 


H a 


Spectral 


Distance 


Adopted 




Photometry 








mag 


mag 


Type 




Luminosity (Lq) 




Preferences 


RY Tau 


04 21 57.40 


+28 26 35.5 


10.2 


6.1 


F8IIIc (1) 


142+14 (2) 


5+2 (3) 


4, 


5, 6, 7 


MWC 480, HD 31648 


04 58 46.26 


+29 50 37.1 


7.7 


6.3 


A5V (1) 


lSlt 2 ^ (8,24) 


15+5 (3) 


5, 


6, 9, 10 


AB Aur 


04 59 41.53 


+40 50 09.7 


7.1 


5.1 


Al (11) 


144±^ (12) 


70+20 (3) 


4, 


5, 6, 7, 9, 10, 13, 1. 


HD 45677 b , MWC 142 


06 28 17.42 


-13 03 11.0 


8.1 


6.3 


B2 (17) 


1000+500 (17) 


14000+7000 (3) 


4, 


18 


MWC 166, HD 53367 


07 04 25.52 


-10 27 15.7 


7.0 


6.2 


BOIVe (19) 


1150 (20) 


100000+50000 (3) 


4, 


5, 15, 18 


HD 144432, HIP 78943 


16 06 57.95 


-27 43 09.4 


8.2 


6.5 


A9IVcv (1) 


145 (21) 


15+4.0 (3) 


4, 


5, 15, 22, 23 


MWC 863, HD 150193 


16 40 17.92 


-23 53 45.2 


8.9 


6.2 


A2IVc (1) 


150±3q (8,24) 


20+10 (3) 


4, 


5, 6, 9, 14 


MWC 275, HD 163296 


17 56 21.29 


-21 57 21.8 


6.9 


5.5 


AlVcpv (1) 


(8,24) 


40+8 (3) 


4, 


6, 9, 16, 22 


MWC 297, NZ Scr 


18 27 39.53 


-03 49 52.0 


12.3 


4.4 


BlVc (1) 


250 (25) 


33000+13000 (3) 


4, 


10, 14, 22 


MWC 614, HD 179218 


19 11 11.24 


+ 15 47 15.6 


7.2 


6.6 


AOV (1) 


240±™ (8,24) 


100+35 (3) 


4, 


5, 6, 15, 22, 26 


vl295 Aql, HD 190073 


20 03 02.51 


+05 44 16.7 


7.8 


6.6 


A2IVcv (1) 


>290 (8,24) 


>83 (3) 


4, 


5, 6, 10 


MWC 342, vl972 Cyg 


20 23 03.61 


+39 29 49.9 


10.6 


5.8 


Bl (27) 


1000 (27) 


33000+15000 (3) 


4, 


5, 27, 28, 29 


MWC 361, HD 200775 


21 01 36.91 


+68 09 47.7 


7.4 


5.5 


B3 (11) + B? (30) 


430±g™ (8,24) 


6000+2000 (3) 


4, 


5, 6, 14 


MWC 1080, v628 Cas 


23 17 25.59 


+60 50 43.6 


11.6 


6.0 


Bl (31) 


2200 (32) 


49000+21000 (3) 


4, 


5, 6, 10, 14 



a Many of the targets are variable stars and these magnitudes (V band from Simbad, and H band from 2MASS) are merely representative. 
b Peculiar B[e] star with uncertain classification 



Note. — References: (1) Mora ct al. (2001), (2) Wichmann et al. (1998), (3) SED fitting, this work, (4) 2MASS; Cutri ct al. (2003), (5) Tycho-2; 
Hog ct al. (2000), (6) USNO-B Catalog; Monet ct al. (2003), (7) Skrutskic et al. (1996), (8) van den Anckcr et al. (1998), (9) Malfait et al. (1998), 
(10) Eisner ct al. (2004), (11) Hernandez ct al. (2004), (12) Bertout ct al. (1999), (13) Herbig & Bell (1988), (14) Hillenbrand et al. (1992), (15) 
Guide Star Catalog 2.2; Space Telescope Science Institute Sz Osscrvatorio Astronomico di Torino (2001), (16) Ducati (2002), (17) dc Winter Sz 
van den Anckcr (1997), (18) New MDM Observatory photometry from 2004 November, this work, (19) Finkenzeller (1985), (20) Herbst ct al. 
(1982), (21) Perez ct al. (2004), (22) DENIS Database; The (2005), (23) Sylvester et al. (1996), (24) ESA (1997), (25) Drew et al. (1997), (26) 
Miroshnichenko ct al. (1999b), (27) Miroshnichenko ct al. (1999a), (28) Kharchenko (2001), (29) Bergner et al. (1995), (30) Pogodin et al. (2004), 
(31) Cohen & Kuhi (1979), (32) Grankin ct al. (1992) 
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Table 2. Observing Log 



Target 


Date 


Interferometer 


Calibrator Names 




(UT) 


Configuration 




RY Tau 


2004 Dec 14 


A28-B10-C00 


HD 27638, HD 27159, HD 26553 


MWC 480 


2002 Dec 13 


A35-B15-C15 


HR 1626 


AB Aur 


2002 Dec 8-14 


A35-B15-C15 


HR 1626, SAO 57504 




2002 Dec 16,19 


A35-B05-C10 


HR 1626 




2003 Feb 19 








zUUo rCD ZU,ZZ,Zo 


A QK TO 1 r ni fi 


rirt lOZO 




onnQ i\T™r on oo 
zUUo J.NOV ZU,ZZ 


Aoo-jdIo-L/IU 


rlrv, 10Z0 




9004 TV<^ 0*3 
ZUU^ L/CC UO 


A ^ R1 ^ POO 


TTR 1 R9R 
nrt ±OZO 




9004 Dpp 1 ^ 


A9K RIO POO 


TIP? 1 R9fi 
ij.iv, j_UiiU 


xlJJ 4DD I ( 


onnQ ivT^,r on on 
zUUo 1NOV Zy,oU 


A QC T3 1 C pin 


rrn /I fiOl Q 
xlJJ 40Zlo 




onn/i A/T^*. ofi 
ZUU4 IVlar ZD 




TTTA /I fiOl Q 

xlJJ 40zlo 




onn/i "Pi^^ 1 q 1/1 
ZUU4 DCC lo, 14 


A oq "Din rvin 
AZo-Jt5iU-L/UU 


un /i sol q 




ZUU4 JJCC Id 


a oq "D n k n 
AZo-JtSUo-L/lU 


TTTTpi A f?Q1 Q 

riiJ 40Zlo 


TV KWT(~ f 1 CtG. 
MWVj 100 


onno int™, on on 
zUUo JNov Zy-oU 


a or dic ni n 
AoO-olD-OlU 


TUT\ cnnC7 TJTA /Id OQ 

riD oUUtw, riJJ 401zo 




9004 Ttr^r- 1 
ZUU4: JJCC 1U 


A 3^ "R1 ^ POO 


un 470^4 TTFl 4Q147 


TTT~"1 1 /I /I /I on 

riD 1444oz 


onn/i T„», o 
zUU4 Jun D-o 


A or dic ni n 


riD 14o/00 


MWC 863 


2003 Jun 17,18 


A35-B15-C10 


HR 6153 




2004 Mar 16,23-25 


A35-B15-C10 


HR 6153 


MWC 275 


2003 Jun 12,21 


A35-B15-C10 


HR 6704, HD 171236 




2004 Mar 26 




HR 6704, HD 171236 




2004 Jun 2,7-8 




HD 157546, HD 170657, HD 174596 


MWC 297 


2003 Jun 12, 23 


A35-B15-C10 


HD 171236, HR 7149, HR 7066 




2005 Jun 20, 22, 26 


A35-B15-C10 


HR 7066, HR 7149 




2005 Jun 29, 30 


A35-B15-C00 


HR 7066, HD 169268 


MWC 614 


2003 Jun 16-17 


A35-B15-C10 


HD 182101, HD 184502 


vl295 Aql 


2003 Jun 14,17 


A35-B15-C15 


HD 185209 


MWC 342 


2003 Jun 13 


A35-B15-C10 


HD 191589, HD 199547 


MWC 361-A 


2003 Jun 21,22 


A35-B15-C10 


HR 7967 




2003 Nov 27-29 


A35-B15-C10 


HR 7967, HD 197950 




2004 Jun 6 


A35-B15-C10 


HD 193664 




2004 Dec 10 


A35-B15-C00 


HR 7967 




2005 Jun 17,18 


A35-B15-C10 


HD 193664 




2005 Jun 29 


A35-B15-C00 


HD 193664 


MWC 1080 


2003 Nov 30 


A35-B15-C10 


HD 221639 




2004 Dec 10 


A35-B15-C00 


HR 8881 




2004 Dec 13 


A28-B10-C00 


HR 8881, HD 221639 



Configuration refers to the location of telescopes A,B,C on the NE, SE and NE arm respectively; see §2.1 for more details 
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Table 3. Calibrator Information 



Calibrator 


Spectral 


Adopted Uniform Disk 


Reference (s) 


Name 


Type 


Diameter (mas) 




HD 26553 


A4III 


0.37±0.30 


getCal a 


HD 27159 


K1III 


0.30±0.10 


gctCal 


HD 27638 


B9V 


0.32±0.10 


gctCal 


HD 46218 


A5 


0.36±0.60 


getCal 


HD 47054 


B8 


0.10±0.15 


gctCal 


HD 49147 


B9.5 


0.35±.03 


gctCal 


HD 50067 


K4III 


1.35±0.35 


gctCal 


HD 143766 


F7V 


0.30 ±0.05 


gctCal 


HD 157546 


B8V 


0.22±0.06 


getCal 


HD 168415 


K4III 


2.6±0.3 


Blackbody fit 


HD 169268 


F6III 


0.28 ± 0.03 


getCal 


HD 170657 


K1V 


0.58±0.10 


gctCal 


HD 171236 


K1III 


1.0±0.5 


gctCal 


HD 174596 


A3V 


0.2±0.05 


gctCal 


HD 182101 


F6V 


0.4±0.1 


getCal 


HD 184502 


B3III 


0.19±0.10 


getCal 


HD 185209 


K3III 


0.7±0.2 


Blackbody fit 


HD 191589 


K5III 


1.47±0.37 


getCal 


HD 193664 


GOV 


0.22±0.07 


gctCal 


HD 197950 


A8V 


0.32±0.06 


gctCal 


HD 199547 


K0III 


0.75±0.15 


getCal 


HD 221639 


K1V 


0.51±0.10 


getCal 


HR 1626 


K0III 


1.2±0.7 


gctCal 


HR 6704 


K0III 


1.89±1.47 


getCal 


HR 7066 


KOI 


2.17±1.01 


gctCal 


HR 7149 


K2III 


1.69±0.75 


gctCal 


HR 7967 


G8III 


0.76±0.30 


getCal 


HR 8881 


M2III 


2.4±0.5 


van Belle (1999) 


SAO 57504 


K0 


0.58±0.36 


gctCal 



''getCal is maintained and distributed by the Michclson Science Center 
(http: / / msc.caltech.edu) 
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Table 4. Herbig Ae/Be Disk Properties 



Target 


Dust Fraction 


Ring Diameter 




Name 


at H-band a 


mas 


AU 


Comments 


RY Tau 


0.8±0.05 


3.99±0.35 


0.57±0.05 


Poor fit w/ no halo emission 






2.85±0.24 


0.40±0.04 


Good fit w/ 10% halo emission 


MWC 480 


0.54±0.06 


3.8±0.3 


0.50±0.04 




AB Aur 


0.65±0.10 


3.2±0.3 


0.46±0.05 


7.5±2.5 % Halo emission 


HD 45677 


0.46±0.03 b 


~15.6 


15.6 


Strong Asymmetry. See Table 6 




0.62±0.03 b 


~12.9 C 


12.9 


Gaussian FWHM 


MWC 166 


O.liO.l 


10±5 


11.5±5.8 


Weak IR excess. Nearby companion star. 


HD 144432 


0.5±0.1 


4.1±1.0 


0.57±0.14 




MWC 863 


0.65±0.15 


4 y+0.8 
4 - '-0.5 


71+ 012 

u - (i -0.08 


Poor SED fit 


MWC 275 


0.59±0.06 


3.3±0.5 


0.40±0.06 


5.0±2.5 % Halo emission 


MWC 297 


0.60±0.05 


5.6±0.5 


1.40±0.13 


Possible <5% Halo 


MWC 614 


0.39±0.03 d 


>13 


>3.1 


Suspect nearby companion 


vl295 Aql 


0.5±0.1 


3.48±0.40 


1.01±0.12 


assuming d = 290pc 


MWC 342 


0.75±0.10 


3.1±0.3 


3.1±0.3 




MWC 361-A 


0.6±0.1 


3.47±0.44 


1.49±0.19 


Primary component only 


MWC 1080 


0.90±0.10 


3.11±0.21 


6.84±0.46 





a Best estimate for fraction of H-band light coming from circumstellar material based on most recent 
photometry. Upper and lower limits are based on SED fitting to diverse data sets and represent the range 
of possible values given historical variability. 

b H-band excess derived directly from fitting visibility data. SED analysis for HD 45677 yields estimate 
for H-band dust fraction 0.65±0.15. 

c Full- Width at Half-Maximum (FWHM) result from fitting Gaussian profile instead of ring model. 

d H-band excess derived directly from fitting visibility data. SED analysis for MWC 614 yields estimate 
for H-band dust fraction 0.37±0.10. 
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Table 5. Closure Phase Results 



Target 


Maximum Closure 
Phase a (deg) 


Minimum Fringe 
Spacing (mas) 


Ring Diameter 
(fringes) 


RY Tau 


0.1±0.4 


11.5 


0.35 


MWC 480 


0.1±2.5 


8.9 


0.43 


AB Aur 


-4.1±0.4 


11.7 


0.27 


HD 45677 


-26.6±4.0 


11.4 


1.37 


MWC 166 


1.8±0.9 


12.1 


0.83 


HD 144432 


-1.6±1.3 


13.5 


0.30 


MWC 863 


5.1±2.6 


11.6 


0.41 


MWC 275 


-0.6±0.4 


11.8 


0.28 


MWC 297 


-1.9±0.4 


9.6 


0.58 


MWC 614 


4.1±1.7 


9.3 


>1.4 


vl295 Aql 


1.9±1.7 


9.4 


0.37 


MWC 342 


-1.8±1.6 


8.9 


0.35 


MWC 361-A 


31. Oil. 5 


10.2 


binary 


MWC 1080 


-1.3±0.5 


9.8 


0.32 



a This column contains the most statistically-significant closure phase datum showing a 
deviation from zero, based on uv-averaged data. The value is bolded if >2-a detection. 



Table 6. Elongated Ring Models for HD 45677 



Model 


Dust Fraction Star Fraction 


Ring Diameter 


Skew 


Skew-PA 


Type 


at H-band 


Major Minor 


PA 




Elliptical Ring 


0.46±0.01 0.54±0.01 


17.8±1.3 12.9±1.4 


70±16 




Skewed Elliptical Ring 


0.46±0.01 0.54±0.01 


18.5±1.5 13.4±0.6 


77±13 0.92±0.08 


-31±8 
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Table 7. Astrometric Binary Solutions for MWC 361- A 



Epoch 


Act 


A5 


Separation 


PA 


(U.T.) 


(mas) 


(mas) 


(mas) 


(°) 


1998 Jun 13-21 


0.9±3.0 


17.9±2.6 


17.9 


2.7 


1998 Sep 29 


2.9±3.1 


17.7±1.1 


18.0 


9.4 


2003 Jun 21-22 


11.4±4.6 


3.5±2.9 


11.9 


73.0 


2003 Nov 27-29 


5.8±3.8 


-5.8±1.0 


8.2 


134.8 


2004 Jun 6 


1.7±2.6 


-11.3±1.2 


11.4 


171.6 


2004 Dec 10 


-3.3±3.4 


-1.9±0.6 


3.9 


239.8 


2005 Jun 17-29 


-2.6±2.0 


12.5±0.8 


12.8 


348.3 



Other binary information: 

H band ratio: 6.5±0.5, Primary UD diameter: 3.6±0.5 mas 
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Table 8. Preliminary Orbital Parameters for MWC 361-A 



Parameter 


Radial Velocity 


Interferometry 




Pogodin et al. 2004 


This Work 


Period (days) 


1341±41 


1377±25 


T (JD) 


2449149±87 


2449152±90 


ai sin(i) (AU) 


1.33±0.08 




a (mas) 




15.14±0.70 


e 


0.29±0.07 


0.30±0.06 


i(°) 




65±8 




203±22 


224±16 


«(°) 




-0.2±7.6 


K (km/s) 


11.2±0.7 




f(M 2 ) (M ) 


0.175±0.035 




M 1 + M 2 (M ) a 


10.4lf 9 5 



a For distance d= 430+gQ pc; see description of system 
mass estimate in §6.2. 
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Fig. 4. — These panels show the observed IOTA3 Closure Phase results as a function of hour 
angle for HD 45677 and MWC 361-A. Each observing configuration and epoch are denoted 
by a unique plotting symbol according to the legend in the left-hand panel. 
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Wavelength (microns) Baseline (m) 



Fig. 5. — (a) left-panel. Spectral Energy Distribution for RY Tau including a simple two- 
component fit. The dashed line represents the stellar contribution only and the solid line 
includes blackbody dust emission. Photometry references included in Table 1. No error bars 
are shown because intrinsic variability dominates over measurements uncertainty at any 
given epoch, (b) right-panel. Visibility data are presented along with best-fit ring models 
exploring different values for halo and dust fractions (see inset legend for details). 
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Fig. 6. — (a) left-panel. Spectral Energy Distribution for MWC 480 including a simple 
two- component fit. The dashed line represents the stellar contribution only and the solid 
line includes blackbody dust emission. Photometry references included in Table 1. No error 
bars are shown because intrinsic variability dominates over measurements uncertainty at any 
given epoch, (b) right-panel. Multiple estimates for the dust fraction are used to constrain 
ring model fits to the IOTA3 visibility data. 
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Fig. 7. — (a) left-panel. Spectral Energy Distribution for AB Aur including a simple two- 
component fit. The dashed line represents the stellar contribution and the solid include 
blackbody dust emission. Photometry references included in Table 1. No error bars are 
shown because intrinsic variability dominates over measurements uncertainty at any given 
epoch, (b) right-panel. Multiple estimates for the dust fraction are used to constrain ring 
model fits to the IOTA3 visibility data; here we also include 7.5% emission from an extended 
halo to fit short baseline data. 
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Fig. 8. — (a) left-panel. Spectral Energy Distribution for HD 45677 including a simple 
two- component fit. The dashed line represents the stellar contribution and the solid include 
blackbody dust emission. Photometry references included in Table 1. Because this SED has 
changed significantly for this target in recent years, we only include most recent photometry 
(see Table 1 for references), (b) right-panel. Here we fitted the IOTA3 visibility with a 
symmetric ring model as well as a Gaussian emission model. Large residuals here are a sign 
of elongated asymmetric emission and a more complicated model is presented in Figure 24. 
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Fig. 9. — (a) left-panel. Spectral Energy Distribution for MWC 166 including a simple 
two-component fit. The dashed line represents the stellar contribution and the solid include 
blackbody dust emission. Photometry references included in Table 1. No error bars are 
shown because intrinsic variability dominates over measurements uncertainty at any given 
epoch, (b) right-panel. Because of a known, nearby companion which is completely resolved 
in our observations, the visibility does not go to unity at the origin. Multiple estimates for 
the companion contribution and dust fraction are used to constrain ring model fits to the 
IOTA3 visibility data. 
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Fig. 10. — (a) left-panel. Spectral Energy Distribution for HD 144432 including a simple 
two-component fit. The dashed line represents the stellar contribution and the solid include 
blackbody dust emission. Photometry references included in Table 1. No error bars are 
shown because intrinsic variability dominates over measurements uncertainty at any given 
epoch, (b) right-panel. Multiple estimates for the dust fraction are used to constrain ring 
model fits to the IOTA3 visibility data. 
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Fig. 11. — (a) left-panel. Spectral Energy Distribution for MWC 863 including a simple 
two-component fit, after removing flux from the T Tauri companion (Corporon 1998). The 
dashed line represents the stellar contribution only and the solid line includes blackbody 
dust emission. Photometry references included in Table 1. No error bars are shown be- 
cause intrinsic variability dominates over measurements uncertainty at any given epoch, (b) 
right-panel. Visibility data are presented along with best-fit ring models covering a large 
uncertainty range in fraction of H band emission arising from disk. 



- 52 - 




Wavelength (microns) Baseline (m) 



Fig. 12. — (a) left-panel. Spectral Energy Distribution for MWC275 including a simple 
two- component fit. The dashed line represents the stellar contribution and the solid include 
blackbody dust emission. Photometry references included in Table 1. No error bars are 
shown because intrinsic variability dominates over measurements uncertainty at any given 
epoch, (b) right-panel. Multiple estimates for the dust fraction are used to constrain ring 
model fits to the IOTA3 visibility data; here we include 5% emission from extended halo to 
fit shortest baselines. 
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Fig. 13. — (a) left-panel. Spectral Energy Distribution for MWC297 including a simple two- 
component fit. The dashed line represents the stellar contribution only and the solid include 
both star and blackbody dust emission. Photometry references included in Table 1. No error 
bars are shown because intrinsic variability dominates over measurements uncertainty at any 
given epoch, (b) right-panel. Multiple estimates for the dust fraction are used to constrain 
ring model fits to the IOTA3 visibility data. 
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Fig. 14. — (a) left-panel. Spectral Energy Distribution for MWC 614 including a simple 
two- component fit. The dashed line represents the stellar contribution only and the solid 
line includes blackbody dust emission. Photometry references included in Table 1. No error 
bars are shown because intrinsic variability dominates over measurements uncertainty at 
any given epoch, (b) right-panel. Visibility data are presented along with best-fit ring and 
Gaussian models. 
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Fig. 15. — (a) left-panel. Spectral Energy Distribution for vl295 Aql including a simple 
two-component fit, after removing flux from the T Tauri companion (Corporon 1998). The 
dashed line represents the stellar contribution only and the solid line includes blackbody 
dust emission. Photometry references included in Table 1. No error bars are shown be- 
cause intrinsic variability dominates over measurements uncertainty at any given epoch, (b) 
right-panel. Visibility data are presented along with best-fit ring models covering a large 
uncertainty range in fraction of H band emission arising from disk. 
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Fig. 16. — (a) left-panel. Spectral Energy Distribution for MWC 342 including a simple 
two- component fit. The dashed line represents the stellar contribution only and the solid 
line includes blackbody dust emission. Photometry references included in Table 1. No error 
bars are shown because intrinsic variability dominates over measurements uncertainty at any 
given epoch, (b) right-panel. Multiple estimates for the dust fraction are used to constrain 
ring model fits to the IOTA3 visibility data. The poor fit might indicate elongated emission, 
however our limited dataset does not allow further investigation at this time. 
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Fig. 17. — (a) left-panel. Spectral Energy Distribution for MWC 361- A including a simple 
two-component fit. The dashed line represents the stellar contribution only and the solid 
line includes blackbody dust emission. Photometry references included in Table 1. No error 
bars are shown because intrinsic variability dominates over measurements uncertainty at any 
given epoch, (b) right-panel. Here, we present the corrected visibility data of the primary 
component after removing the close binary companion. No formal error bars are shown 
since the extraction process introduces an uncertain error, however the scatter in values 
should represent the uncertainty. Visibility fits are presented along with best-fit ring models 
covering a large uncertainty range in fraction of H band emission arising from disk (based 
on SED decomposition). 
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Fig. 18. — (a) left-panel. Spectral Energy Distribution for MWC 1080 including a simple 
two- component fit. The dashed line represents the stellar contribution only and the solid 
line includes blackbody dust emission. Photometry references included in Table 1. No error 
bars are shown because intrinsic variability dominates over measurements uncertainty at any 
given epoch, (b) right-panel. Visibility data are presented along with best-fit ring models 
covering a large uncertainty range in fraction of H band emission arising from disk. 
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Fig. 19. — This figure shows the telescope geometry for closure phase measurement. Turbu- 
lent air pockets cause optical path length variations above individual telescopes corrupting 
fringe phase measurements. However by summing fringe phases around a closed triangle, 
these phase errors cancel out resulting in a good interferometric observable, the closure phase. 
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Fig. 20. — This figure shows an example closure phase calculation for a DDN disk model 
using IOTA, (a) The top left panel shows the near-IR synthetic image assuming an inner 
rim radius of 3 milliarcseconds, viewed at inclination % = 30°. (b) The top right panel shows 
a typical observing configuration for the IOTA 3-telescope interferometer, (c) The bottom 
left panel shows the closure phase as a function of hour angle, assuming object at declination 
5 = 35° (dashed line for the DDN model rotated 90° on the sky), (d). The bottom right 
panel shows the closure phase signal for baselines twice as long, illustrating sensitivity with 
angular resolution. 
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Fig. 21. — This figure shows predicted IOTA closure phase signal for examples of "skewed 
ring" models with similar geometry to the DDN model seen in Figure 20 (see §4.4.2 for model 
details), (a) Top-left panel shows synthetic image for ring model with maximum skew = 1. 
(b) Top-right panel shows synthetic image for ring model with intermediate skew = 0.5. (c) 
Bottom-left panel shows the closure phase result for skew = 1 disk (same IOTA observing 
geometry as in Figure 20b). (d) Bottom-right panel corresponds to ring model with skew 
= 0.5. 
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Fig. 22. — In this figure we compare the maximum predicted closure phases for DDN models 
as a function of disk inclination angle and angular resolution of interferometer. For well- 
resolved disks, the closure phase signals for DDN models are expected to be huge but are 
not observed in any targets of our sample. A few of the most notable targets are labeled. 
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Fig. 23. — In this figure we compare the maximum predicted closure phases for skewed 
ring models as a function of the skew amount and angular resolution of interferometer. 
Note that the line for skew =0 is identically zero due to centrosymmetry. Note that most 
systems are not resolved enough to constrain the skewed models, calling for longer-baseline 
interferometers such as PTI, CHARA, and VLTI. 
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Fig. 24. — Top-left panel shows the best-fit Elliptical Ring Model to the IOTA3 data (North is 
up, East is left). Top-middle panel compares the observed visibility data with the model data. 
Top-right panel compares the observed closure phases with the model closure phases. Since 
this model is constrained to be centrosymmetric, all model closure phases are identically zero 
and thus a good fit is impossible. The bottom panels present the best fit "skewed elliptical 
ring model." The closure phases are fairly well-explained by a strong north-west skew. All 
parameters of these fits are presented in Table 6. Some data points are not well-explained 
by these models suggesting a patchier dust distribution. 
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Fig. 25. — This figure shows the astrometric orbit for MWC 361-A (solid line). The observed 
locations of the secondary component (relative to primary) are marked with error ellipses 
(dotted lines with crosses) for seven different epochs. The dashed line connects the locations 
of the primary and secondary star at periastron passage. 



